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Abstract. Boolean networks have been successfully used in modelling gene regu- 
latory networks. In this paper we propose a reduction method that reduces the 
complexity of a Boolean network but keeps dynamical properties and topological 
features and hence it makes the analysis easier; as a result, it allows for a better 
understanding of the role of network topology on the dynamics. In particular, we 
use the reduction method to study steady states of Boolean models. 



1. Introduction 

Boolean networks have been successfully used in modelling gene regulatory networks 
such as the Drosophila segment polarity network [l|, the yeast cell-cycle network Q| and 
the Th regulatory network [jj. Boolean networks provide a nice theoretical framework 
that allows for simulation, control theory and reverse engineering. 

However, their analysis is not a trivial task. For example, the problem of finding 
steady states has been shown to be NP-complete [llj. One way to overcome these kind 
of problems is to develop mathematical and computational tools 0, 0> EL @ • While these 
tools allow to answer some questions, such as what the steady states are, it is often not 
intuitive why such answers were obtained. Another way is to reduce the network to one 
that has less complexity while keeping the main features; the reduced network is easier to 
analyze and can not only help to answer questions, but also to give insight of why such 
answers were obtained. This in turn, provides a better understanding of the problem that 
is being studied. 

In this paper a reduction method for Boolean networks is proposed; a preliminary 
formulation of this method was provided in fici| ]. The reduction method reduces the 
complexity of Boolean networks making the analysis easier and also elucidates the role of 
network topology in dynamics. We will focus on the existence, number and type of steady 
states. A similar reduction method for logical models has been proposed in [7|. 

This paper is organized as follows, in Section [2] we present the reduction method; 
properties are presented in Section [3] Section [4] includes an application of the reduction 
method. We close with a discussion in Section [S] 



2. Reduction Method 

2.1. Reduction Steps. We now provide the reduction steps to reduce a Boolean net- 
work and its corresponding wiring diagram. The idea behind the reduction method is 
simple: the wiring diagram and Boolean functions should reflect direct regulation and 
hence nonfunctional edges and variables should be removed; on the other hand, vertices 
can be deleted, without losing important information, by allowing its functionality to be 
"inherited" to other variables. 

(1) We simplify the Boolean functions and wiring diagram: 

(a) Reduce Boolean expressions using Boolean algebra. This will delete variables 
that are not functional. 

(b) Delete edges that do not correspond to Boolean expressions. That is, we 
delete edges that are non functional. 
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(2) We delete vertices with no self loop, that is, vertices whose Boolean function does 
not depend on it. Let x% be a vertex such that f Xi does not depend on Xi. 

(a) For all vertices Xi — > y, that is, for all vertices whose Boolean function 
depends on Xi, replace the Boolean function for y, f y (xi, . . . , Xi, . . . , Xh), by 

fy{X\ , • ■ • , fxi ) ■ • ■ 3 ^fc)- 

(b) Replace edges y^Xi-^z by y~*z and delete Xi (and edges from/to Xi) 
We will see that these steps give rise to a reduced network that keeps features of the 

wiring diagram and dynamical properties of the original network. 

2.2. Reduction Algorithm. We present an algorithm to simplify Boolean functions and 
their wiring diagram (S); and an algorithm to eliminate a vertex x (R) . 

2.2.1. Algorithm S. Input: / = (/i, . . . , /„) and A. 

(1) For i = 1, . . . , n: 

(2) Simplify fi using Boolean algebra 

(3) Construct A corresponding to variables appearing in fi 
Output: (Simplified) / = (/i, . . . , /„) and A. 

Example 2.1. Consider the Boolean network given by / = (/i>/2,/3) = ((x2 A X3) V 
X2, (xi A X3) V -1x2, -UEi) with wiring diagram given in Figure[T] Algorithm S gives as an 
output / = {fi,fz,fs) = (x2, (xi A13) V -1X2, ->xi) with wiring diagram given in Figure[T] 
We can clearly see that Algorithm S detected that the although the variable x$ appears in 
fi , it is not functional; hence X3 is removed from /1 (using Boolean algebra) and the edge 
xz — > xi is deleted as well. That is, after using S, we obtain an accurate representation of 
/ and its wiring diagram. 




Figure 1. Wiring diagram of / before (left) and after (right) using 
algorithm S. Arrows indicate positive paths and circles indicate negative 
paths. 

If we say that a variable z depends on x, we mean that f z depends on x. Also, if we 
write fi = fi{xi ± , . . . , Xi k ), we are saying that we are considering the function fi in terms 
of Xi 1 , . . . , Xi k , even if fi does not depend on some of the Xi . 's; however, it does mean 
that fi does not depend on the other variables. This convention will make some of the 
definitions and proofs simpler. 

Let / = (/1, . . . , /„) be a Boolean network with wiring diagram given by the adjacency 
matrix, A. We show the algorithm to reduce / by eliminating a vertex x that does not 
have a self loop. The algorithm to eliminate x is as follows: 

2.2.2. Algorithm R. Input: / = (/1, . . . , /„), A and x that does not depend on itself. 

(1) Find the variables that depend on 1: zi, ■ ■ ■ ,z r 

(2) Fori = l,2,...,...,r 

Replace f Zi = f Zi (x,j . .) by f H = f H (f x ,.. .) 

(3) Let f - x - — (/1, . . . , f x , . . . , /„) (where f x means that f x is omitted) and simplify 
it using S. 

(4) Let A^'^adjacency matrix of the wiring diagram of f^ 
Output: f [x] and A [x] . 
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Example 2.2. Consider the Boolean network given by / = (/i, fa, fa) = (2:2, (xi A X3) V 
-1X2,-1X1) with wiring diagram given in Figure [2] After using algorithm R we obtain the 
network j' 13 ' = h = (fti, fo) = (2:2, (xi A -1x1) V -13:2) = (X2, ~^x 2 ) with wiring diagram in 
Figure[2]. We can see that the functionality of £3, that is, "being -1x1", is inherited to x 2 
and simplified using S. 




Figure 2. Wiring diagram of / before (left) and after (right) using 
algorithm R to eliminate 2:3. 

The Boolean network obtained by eliminating vertices Xi x , . . . , Xi k is denoted by f[ x *i>---> Xi k 
with wiring diagram A^ Xi i '■■■> Xi k\ The following proposition states that the order in which 
variables are eliminated in not important. 

Proposition 2.3. Using the notation of Algorithm R we have the following: Using the 
notation above we have the following: f^ x " x i^ = f\*j> x d and A lxi,x ^ = A lx j' Xi ^ 

Proof. Without loss of generality we can assume Xi = Xi and Xj = x n . 

Let g' = (g[, . . .,9'^) = /W, g = (g 2 , . . .,g n -i) = g'^ = /t»».»il, h' = (h' 2 , ...,ti n ) = 
f lxi] , h= (h 2 ,...,h n -i) = h! [xn] =/[«!.*»]. We consider 4 cases. 

Case 1. x n and x\ do not depend on each other. 

Then, since fi does not depend on x n (nor itself), we have /1 = fi(x 2 , ■ ■ ■ ,x n -i) and 
since /„ does not depend on xi (nor itself) we have /„ = f n {x 2 , ■ ■ ■ ,x n -i). 

Then, g[ = /i(2;i,2;2, . . . , x n _i) and g\ = fi(xi,x 2 , ■ ■ ■ ,x n -i, f„(x 2 , . . . ,x n _i)) for i = 
2, . . . , n — 1. Then, for i = 2, . . . , n — 1 we have 

9i = g'i{g'i,x 2 , . . .,X n -l) 

= fi(9l,X2, ■ ■ ■,X„-l,f„(x 2 , ■ ■ .,X n -l)) 

= fi(fl( x '2, ■ • ■ , X n -\),X 2 , . . . ,X n -!, f„(x 2 , . . . , X n -i)) 

Similarly, it follows that hi = fi(fi(x 2 , ■ ■ ■ , x n -i), x 2 , . . ■ , x n -i, fn(x 2 , . . . , a;„_i)) for 
i = 2, . . . , n — 1. Hence, gt — hi for i = 2, . . . , n — 1 and g = h. 

Case 2. x n depends on x\ but xi does not depend on x n . 

Then, since /„ depends on xi (and not on x n ), we have /„ = f n (xi, . . . ,x n -i) and 
since fi does not depend on x n (nor xi), we have /1 = fi(x 2 , . . . , x n -i). 

Then, g[ = fi(x 2 , . . . ,a;„-i) and^ = fi(xi, . . . , x n -i, fn(xi, . . . ,x n -i)) for i = 2, . . . , n- 
1. Then, for i = 2, . . . , n — 1 we have 

9i = g'i(gi,x 2 , ■ ■ -,x n -i) 

= fi(g'i,x 2 , ■ . . ,x n - 1 ,f n (g' 1 ,x 2 , . . .,x n -i)) 

= fi{h(x2, ■ ■ -,X n -l),X2, ■ ■ . ,X n -i, f n (f!(x 2 , . . .,X„-l),X2, ■ ■ .,X n -i)) 
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On the other hand, h'i = fi(fi(x2, ■ ■ ■ ,x n -i),X2, 
fn{fi{x2, ■ ■ .,x n -i),X2, ■ ■ .,x n -i)- Then, for i = 2, 



,x n ) for i — 2, . 
, n — 1 we have 



, n — 1 and ft4 = 



ft; = ft;(cc2, . . • , x n -i,h' n ) 

= fi(h(X2, ■ ■ -,X„-l),X2, ■ ■ .,X„-l,h' n ) 



= fi{fl{X2, ■ ■ ■ ,X n -l),X2, ■ ■ ■ , X„-i, fn{fl{X2, X n -l),X 2 , ■ ■ ■ , X n -l)) 

Hence, g; = ft; for i = 2, . . . , n — 1 and g = ft. 

Case 3. xi depends on x„ but x n does not depend on xi. It is analogous to Case 2. 
Case 4. xi and x n depend on each other. 

Then, since /„ depends on x\ (and not on x n ), we have f n = f n (xi, . . . , x n -i) and 
since fi depends on x n (and not on xi), we have /i = fi(x2, ■ ■ ■ ,x n ). 

Then, g[ = /i(a; 2 , . . . , x n -i, f n (xi, . . .,x„-i)) and g[ = fi(xi,. . . , x„-i, f„{xi, . . . ,x n -i)) 
for i = 2, . . . , n — 1. Then, for i = 2, . . . , n — 1 we have 

5» = Si (fli, • • -,Xn-i) 

= fi(g'l,X2, • • ■ ,£„-!,/„(#!, . . . ,X n _l)) 

= fi(fl( x 2, ■ ■ ■ , X n -1, fn{xi, . . . ,X„-l)),X2, ■ ■ ■ , X n -1 , fn (/l {X2 , ■ ■ ■ , X n - 1 , fn (xi , . . . , X n _l)), X2, . . . , X n - 

On the other hand, h' n = fn(fi(x2, ■ ■ ■ ,x„),X2, x„-i) and h'i = fi(fi(x2, ■ ■ ■ ,x n ),X2, ■ ■ ■ ,x„) 
for i = 2, . . . , n — 1. Then, for i = 2, . . . , n — 1 we have 

ft; = /l; (X2, . . . , X n -1 , h' n ) 

= fi(fl(X2,-- .,h' n ),X2, • • • , X n -1, h' n ) 

= fi(fl( x 2, ■ ■ ■ , Xn-1, fn{fl{X2, . . . , X n ) , X2, ■ ■ ■ ,X n -l)),X2, ■ ■ ■ , X n -1 , fn (/l {X2, ■ ■ ■ , X n ), X2, ■ ■ ■ ,X n -l)) 

Notice that fi(x2, ■ ■ ■ , x n -i, f n (xi, . . . , x n -i)) does not depend on xi and f n (fi {x2, ■ ■ ■ , x n ), X2, ■ ■ ■ , x n -i) 
does not depend on x n (g = g'^ and h = ft'' 1 "' would be undefined otherwise). Then, 
for i = 2, . . . , n — 1 

9i = fi{fl( x 2, ■ ■ ■ , Xn-1, fn(xi, . . . ,X n -l)), X2, ■ ■ ■ , X n -1, fn (/l {X2, ■ ■ ■ ,X n -l, fn(xi, . . . ,X n -l)),X2, ■ ■ ■ ,X n - 
= fi{fl( x 2, ■ ■ ■ , Xn-1, fn{xi, . . . , X n -l)),X2, ■ ■ ■ , X n -1, fn (/l {X2, ■ ■ ■ , X n -1, X n ) , X2, ■ ■ ■ ,X n -l)) 
= fi(fl(x2, ■ ■ ■ ,X n -l,fn(fl{X2, ■ ■ ■ ,X n ), ■ ■ ■ ,X n -l)),X2, ■ ■ ■ , Xn-1, fn(fl{X2, ■ ■ ■ , X n -1, X n ) , X2, ■ ■ ■ ,X n -l) 
= hi 



The Boolean network obtained by eliminating all variables that can be eliminated 
(variables that do not have a self loop) is denoted by f R with wiring diagram A R . From 
Proposition [] , it follows that f R and A R are independent of the order chosen to eliminate 
vertices (but they do depend on the choice of variables to be eliminated). Also, it may be 
the case that f R and A R are empty; in the case they are not empty, each vertex has a self 
loop. 

Example 2.4. Consider the Boolean network / defined by: 



Hence, g = ft 



□ 



h 
h 

h 
h 
h 
h 



xi A -ii-io 
X2 A -iXlQ 
x :i A -iX8 



Xg V Xn 
->X7 A X\2 





h = x 4 
fs = x 5 



fg = X 6 A -i£l2 

fio = x 7 V xn 



fn = {x 7 V xn) A ->xi2 
fl2 = x 8 A -ixn 
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This Boolean network corresponds to the logical model for Th-lymphocyte differentia- 
tion presented in [jj. It turns out that we can eliminate the variables xi, xi, xz, X4,, x^,x^,xi, xg, xg, xio. 
Then f R = h = (hn,hi 2 ) — (xn A -1X12, _l Xn A X12) ■ Notice that each vertex in the the 
wiring diagram of the reduced network has a self loop. 

3. Properties of the Reduction Method 

3.1. Reduction Method and Dynamical Properties. We now show that the original 
and reduced network share important dynamical properties. The next theorem states that 
the reduction method does not create nor destroy steady states. 

Theorem 3.1. Let f be a Boolean network and g = y ? t !c »i' !c *2'---' x *fcl m th k < n. Consider 
the projection n : {0, l} n —* {0, l} n ~ fc defined by 

w(xi, ... , x n ) = (xi, . . . , x il - 1 ,Xi 1 ,Xi 1+1 , . . . ,Xi h -i,Xi k ,Xi k +i ... , x n ). Then, tt defines a 
one to one correspondence between the set of steady states of f and the set of steady states 

of g- 

Proof. We only need to prove the theorem for k = 1; the general case follows by induction. 
Without loss of generality we can assume that g = (gi, . . . ,g n -i) = f^ Xn \ 

If z = (zi, ■ ■ ■ , Z n ) is a steady state of /, that is, f(z) = z, or /;(z) = Zi for i = 1, . . . , n, 
we want to show that n(z) — (zi, ■ ■ ■ , Zn-i) is a steady state of g, that is, g(n(z)) = n(z) 
or gi(n(z)) = Zi for i = 1, . . . , n — 1. On the other hand, if z = (z\, . . . , z n -i) is a steady 
state of g, that is, g%{z) = Zj for i = 1, . . . , n — 1, we want to show that there is a unique 
steady state of /, z, such that n{z) — z. 

Without loss of generality, suppose x n depends on x% for i = l,...,x r ; and Xi for 
i — s, . . . ,n — 1 are the variables that depend on x n . Notice that it may be the case 
that r = or s — n; in that case x n would depend on no variables or no variable 
would depend on x n . Then, /„ = f„(xi, . . . , x r ); also, for i = s, ...,n — 1 we have 
fi = fi(xn,Xi2, ■ ■ ■ ,x„). Then, for i = 1, . . . , s — 1 we have gi(n(x)) = fi(x) and for 
i = s, . . . ,n — 1 we have g t = fi(xn,Xi 2 , ■ ■ ■ , f n {x\, x r )). 

Let z be a steady state of /. Then, for i = 1, . . . , s — 1 we have <?;(7r(z)) = fi(z) = Zi; 
for i = s,...,n- 1, gi(ir(z)) = fi {zn,Za, . . . , f„(z ly . . . , z r )) = fi(zn,Zi2, ■ ■ ■ , z n ) = Zi. 
Then, g(n(z)) = n(z). 

Let Jbea steady state of g. Define z — (z, f n (zi, . . . , z r )), that is, z n = f n (zi, ■ ■ ■ , z r ); 
we claim that z is a steady state of /. Notice that n(z) — z. For i — 1, . . . , s — 1 
we have that fi(z) = gi(ir(z)) = gi(z) — Zi\ also, for i = s, ...,n — 1 we have that 
fi(z) = fi(za,Zi2, ■ ■ ■ , z n ) = fi(zn,Zi2, • • • , fn{zi,. . . , z r )) = gi(Tv(z)) = Zf, also, f„(z) = 
f n (zi, ■ ■ ■ , z r ) = Zn- This shows that z — (z, f n (zi, . . . , z r )) is a steady state of /. We now 
show that z is unique. Let z' = (z[, . . . , z' n ) be another steady state of / such that tt(z') — 
z; it follows that z = (z, z' n ). Since z' n = f„{z') — f n (z[ , . . . , z' r ) = f n (zi, ... ,z r ) = z n , 
then, z = (z, z n ) = z'. Hence, the steady state is unique. □ 

Corollary 3.2. Let f be a Boolean network and g = f R not empty. Then, there is a one 
to one correspondence between the set of steady states of f and the set of steady states of 

g- 

Corollary 3.3. Let f be a Boolean network and g = f R . If g is empty, then f has a 
unique steady state. 

Proof. Suppose g is empty. Without loss of generality, suppose x\ is the last variable 
to be eliminated. That is, h = f^ 2 ' ■■■• Xn ^ and g = h} Xl \ Since the wiring diagram of 
h : {0, 1} — » {0, 1} has only the vertex x\ and it cannot have a self loop, it follows that 
h cannot be the Boolean function h(x\) — x\ or h(x\) = -ubi. Then, h has to be the 
Boolean function h(xi) = or h(xi) = 1 and hence it has a single steady state (xi — or 
xi — 1, respectively). □ 
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Corollary 3.4. Let f be a Boolean network and g = p x n> x >2''"' x >k' with k < n. Then, f 
is a oscillatory system (the only attractors it has are periodic orbits) if and only if g is a 
oscillatory system. 

Proof, f is an oscillatory system if and only if / does not have any steady state if and 
only if g does not have any steady state if and only if g is an oscillatory system. □ 

3.2. Reduction Method and Topological Properties. We now show that the original 
and reduced network share topological properties. The next theorem states that the 
reduction method does not create new paths nor it changes their signs. 

Theorem 3.5. Let f be a Boolean network and g — p x n' Xi 2>--> x i^ with k < n. Then, if 
there is path from y to z in the wiring diagram of g, there is also a path from y to z in the 
wiring diagram of f (or equivalently, if there is no path from y to z in the wiring diagram 
of f, there is no path from y to z in the wiring diagram of g ). Furthermore, if all paths 
involved in the reduction from y to z in the wiring diagram of f are positive (negative), 
the corresponding paths from y to z in the wiring diagram of g are positive (negative). 

Proof. We only need to prove the theorem for k = 1; the general case follows by induction. 
Without loss of generality we can assume that g = (gi, . . . ,g n -i) = f[ Xn >. Suppose there 
is a path from y to z in the wiring diagram of g. Without loss of generality we can suppose 
the path is y = xi — * X2 — ► • • • — ► x r —i — * x r = z. It follows that gi depends on x%-\ for 
i = 2,...,r. 

We claim that there is a path from xi to X2 in the wiring diagram of /. If X2 depends on 
x„, g2 = fi{xj 1 , . . . ,Xj t , /„). Then, since gi depends on xi, it follows that Xi is one of the 
Xj l 's or f n depends on x\\ then, it follows that there is a path from x\ to X2 in the wiring 
diagram of /. If, on the other hand, x 2 does not depend on x n , g2 = f 2 {xj 1 , . . . , Xj t ). 
Then, since 32 depends on xi, it follows that x\ is one of the Xj^s; then, there is a path 
from xi to X2 in the wiring diagram of /. Similarly, it follows that there is a path from 
Xi-i to Xi for i = 2, . . . , r in the wiring diagram of / and hence, there is path from y — x\ 
to z — x r in the wiring diagram of /. 

Now, to conclude the proof of the theorem it is enough to show that if the paths 
x n -2 — * x n -i and x„- 2 —* x n — > x n -± in the wiring diagram of / = (fi,...,f n ) 
are positive, so is the path x\ — > X3 in the wiring diagram of g = (gi, . . . ,p n _i). 
Since f„(x) = /„(..., x n -2) and f„-i(x) = . . , x n - 2 ,x n ), then g n -i(- ■ ■ , x n -i) = 

f n -\{. . . ,x n -2,f n (. ■ ■ ,x n -2))- To show that the path x„_2 — > x n -i in the wiring di- 
agram of g is positive, we need to show that g n _i(...,0) < g n -i(- ■ ■ , 1)- Since the 
paths x n -2 —> x n -i and x n -2 ~ * x n ~+ x n -\ in the wiring diagram of / are positive, 
then /„_!(..., 0, /„(..., 0)) < /„_!(..., 0, /„(..., 1)) < /„_i(..., 1, /„(..., 1)). Hence, 
g n — 1(. . . , 0) < g n -i(- ■ ■ , 1)- Then, the last part of the theorem follows by induction. □ 

Corollary 3.6. Let f be a Boolean network and g — f R . Then, if there is a path from y 
to z in the wiring diagram of g, there is also a path from y to z in the wiring diagram of 
/■ 

Corollary 3.7. Let f be a Boolean network and g = f\ x i 1 > x i 2 '---> x i k } with k < n. Then, if 
there is a feedback loop at y in the wiring diagram of g, there is also a feedback loop at y 
in the wiring diagram of f . 

Corollary 3.8. Let f be a Boolean network and g — f . Then, if there is a self loop at 
y in the wiring diagram of g, there is also a feedback loop at y in the wiring diagram of f. 

Example 3.9. The next example shows how the reduction method can allow to detect 
nonfunctional paths. Consider the Boolean network: 

/ = (/ii /2; /3i fi) = {pXx,X3 V —<x&, Xi, Xi) with wiring diagram shown in Figure[3] If we 
reduce the network by deleting vertices X3, X4 we obtain f^ x ^< x ^ = g = (<?i, 32) = (~<xi, 1). 
We can clearly see that the reduced network does not have any path from xi to X2. Also, 
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the reduced network does not have any steady state; hence, by Theorem 13.11 / does not 
have steady states. 















(x4)— 


(x2) 




[x2] 



Figure 3. Wiring diagram of the full (left) and reduced (right) network 
for Example 13.91 

Example 3.10. The next example shows how the reduction method can allow to detect 
nonfunctional feedback loops. Consider the Boolean network: 

/ — (/i> /zi /3i /-i) = (^2,2:3 V X4, xi, ->Xi) with wiring diagram shown in Figure [3] We 
can see that there are two feedback loops (one positive and one negative). If we reduce 
the network by deleting vertices £3, X4 we obtain J^ 3 '^ 4 ' = g = (51,52) = ( x 2, 1). We can 
clearly see that the reduced network does not have any feedback loop. Also, it is easy to 
see that g has a unique steady state, (1,1); hence, by Theorem 13. II f has a unique steady 
state. 




Figure 4. Wiring diagram of the full (left) and reduced (right) network 
for Example |3~Tu1 

4. Application 

We now consider the Boolean model presented in Q. It is a small model for Th- 
lymphocyte differentiation. Its wiring diagram is given in Figure |S] The variables and 
Boolean functions of the model are given as follows: 



variable 


Boolean function 


:ci=LNF-7 


fl — Xg V £11 


x 2 =IL-4 


h = -^x 7 A K12 


x 3 =IL-12 


/a = 


x 4 =IFN-7 R 


fi = Xl A -UElO 


x 5 =IL-4R 


fs = X 2 A -UB10 


z 6 =IL-12R 


fe = xg A ->xs 


x 7 =STATl 


h = x i 


z 8 =STAT6 


fs = x 5 


z 9 =STAT4 


fg = xe, A -1x12 


kio=SOCS1 


/10 = x 7 V in 


xu=T-bet 


/11 = (x7 V xn) A -1X12 


xi2=GATA-3 


/12 = isA -<xii 



8 



ALAN VELIZ-CUBA 




Figure 5. Wiring diagram for the Th-Lymphocyte differentiation network. 
Notice that this Boolean network has 2 12 = 4096 states. 

In Example 12.41 we observed that the reduction method gives the following reduced 
model with wiring diagram given in Figure [6] 



variable 


Boolean function 


a;ii=T-bet 


hll = £11 A -iXi2 


zi2=GATA-3 


hl2 = -iXll A 2i2 




Figure 6. Wiring diagram for the reduced network of the Th- 
Lymphocyte differentiation network. 

Notice that the reduced network has only 2 2 = 4 states which is about 0.1% of the 
number of states of the full network. It is easy to see that there are 3 steady states: 
(0, 0), (0, 1), (1, 0). Furthermore, we can explain the nature of the steady states: if 
212 =GATA-3=0, then xn =T-bet can be or 1 for (211,212) to be a steady state; on 
the other hand, if 212 =GATA-3=1, then 211 = T-bet must be 0. Since the reduced net- 
work has 3 steady states, by Theorem l3.ll the larger network has also 3 steady states: si = 
(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), s 2 = (0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1), s 3 = (1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0). 
In the reduced network, the existence of the self loops at GATA3, T-bet and the positive 
feedback loop T-bet<=^ GATA3 suggests that they are the key to the dynamical proper- 
ties; this in turn suggests that the corresponding feedback loops in the larger network are 
determining factors in the dynamics. 

Now, let us see how the reduced network can help us understand the larger network. 
For the reduced network, it is not difficult to see that deleting one or both of the loops at 
GATA3 or T-bet results in the loss of the steady state (0,0). On the other hand, deleting 
one of the edges T-bet ^ GATA-3, GATA-3-^ T-bet does not change the steady states; 
however, if we delete both edges, a fourth steady state (1,1) is created. It is important to 
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notice that this information can be easily obtained. We can expect that the larger network 
has similar properties; to check this we study the effect of deleting edges: deleting the 
loop at T-bet and other edges so that we do not have a feedback loop at T-bet, or any 
edge in the feedback loop [IL-4,IL-4R, STAT6,GATA-3] results in the loss of the steady 
state si that corresponds to (T-bet, GATA3) = (0,0). On the other hand, deleting one of 
the edges T-bet— * GATA-3, GATA-3^ T-bet does not change the steady states; however, 
deleting both edges and other edges so that we do not have a path from T-bet to GATA3 
and GATA3 to T-bet, results in the creation of a fourth steady state corresponding to 
(T-bet, GATA3)=(1,1) . All these properties of the larger network are consistent with 
those of the reduced network that only has 0.1% of the number of states. In summary, the 
reduction method generated a small network that allowed to easily study the existence 
and type of steady states and the role of the feedback loops in the dynamics. 

5. Discussion 

Boolean networks have been successfully used in modeling; they provide a theoreti- 
cal framework that allows for simulation, control theory and reverse engineering. Since 
their analysis is not a trivial task many mathematical and computational tools have been 
developed. 

In this paper we have proposed a reduction method that although simple, it can make 
the analysis of Boolean networks easier. In particular, we applied the reduction method 
to analyze the steady states of a Boolean model for the Th-lymphocyte differentiation. 
The reduction method was not only able to make the analysis of steady states easier but 
was also able to explain the role of feedback loops. 

Future work is to study how the reduction method can help in the analysis of limit 
cycles of a Boolean network. Also, another future project is the generalization of the 
reduction method to general finite dynamical systems. 
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